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Abstract. The expected Riesz energies E n£s of the zero sets of systems 
of independent Gaussian random polynomials are determined asymptotically 
as the degree Af — ^ oo in all dimensions and codimensions. The asymptotics 
are proved for sections of any positive line bundle over any compact Kahler 
manifold. 



1. Introduction 

A classical problem is to determine the optimal configurations of N points 
{z^ , . . . , z^} on the unit sphere S"^, i.e. configurations which minimize some energy 
function, 

(1) £:,(^i,... ,z^)=^G,(z\z^). 

Natural choices of the kernel Gs{z, w) are the Riesz kernels, 

{rg{z,wy, < s < 2m 

~logrg{z,w), s = 0. 

Here, rg{z,w) is the geodesic distance between z,w. The assumption < s < 2m 
ensures that the Riesz function rg{z,'w)~'^ is integrable. The minimal energy is 
denoted by 

(3) £s{N)^ min . . . , z^). 

{zi,...,zjv} 

For background on this problem, we refer to |BBGKSl IHS| . The problem of min- 
imizing energy is of interest on any compact Riemannian manifold, although it 
seems mainly to have been studied on the standard sphere. 

Of course, it is a very demanding problem to determine the actual optimizers. 
Related problems are to find the asymtotics of the minimial energy ([3]) as AT -> oo 
and to determine the properties of asymptotically minimizing configurations, which 
are sometimes called elliptic Fekete points. For instance, to minimize energy the 
points must be well-separated and equidistributed with respect to the volume form. 

This article is concerned with the expected energies E^nEs of the zero sets of 
Gaussian random systems of holomorphic polynomials and their generalizations 
to any compact Kahler manifold. Thus, we do not study optimal configurations 
but rather consider the separation and equidistribution properties of configurations 
defined by random zero sets as measured by the energies M. 
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Our results generalize those of a recent pair of articles [Zl lABSj in the complex 
one dimensional case, which showed that the zeros of Gaussian random holomorphic 
polynomials of one complex variable are asymptotic energy minimizing configura- 



tions to a certain order in N. The results are reviewed in ^1.2 Such polynomials 
are holomorphic sections of 0{N) — > CP^ = 5*^, the iVth power of the hyperplane 
bundle over the Riemann sphere. It has been known for some time that in complex 
dimension one, the zeros of random polynomials of degree N (or of any holomorphic 
sections) repel each other so that the nearest neighbors are almost always of order 
apart. More precisely, the probability of pairs being apart can be deter- 
mined from the scaling asymptotics of the two-point correlation function of zeros 
in [Hi IBSZlT IBSZ21 rBSZ3| . Furthermore, the 'empirical measure of zeros' defined 

by 

(4) = ^ E 

z:s«(z)=0 

tends to the 'equilibrium measure' given by the Riemannian area form [SZ2]. The 
leading order asymptotics of the energy mainly reflect this repulsion and equidis- 
tribution. Repulsion and equidistribution are not sufficient to ensure that random 
zero point configurations asymptotically minimize energy, and indeed one would 
not expect a random configuration to do as well as a carefully chosen one. But in 



complex dimension one, random zero sets do surprisingly well (see ^1.2). 

In this article, we consider the expected energies of the zero ses of a random full 
system of m holomorphic polynomials on a complex projective space CP™ of any 
dimension. In fact, we consider the problem for sections of ample holomorphic line 
bundles L M over general Kahler manifolds, and moreover for systems of fc < m 
holomorphic sections. We refer to the case of a full system of m polynomials in 
complex dimension m, whose simultaneous zero set is (almost surely) a discrete set 
of points, as the point case or equivalently as the discrete case. 

In higher dimensions, we do not expect the complex zeros in the discrete case to 
asymptotically minimize energy. They do become equidstributed with respect to the 
equilibrium measure but they are not well separated. As recalled in §2.3[ the scaling 
limits of the pair correlation function of zeros were calcuated in |BSZlllBSZ2[|BSZ3] 
in all dimensions and co-dimensions. The results show that, in the point case, the 
zeros almost ignore each other with only a slight repulsion at very close distances 
in dimension 2, and they in fact attract or clump together in higher dimensions. 
The results of this article give a somewhat different perspective on the separation 
of the zeros as reflected by energy asymptotics, or more precisely on the distance of 
zeros from the discriminant variety of polynomial systems with a double root (see 



^2.41 



Before stating our results on expected energies, let us review what is known 
about the minimum energies of configurations of points in the point case so that 
we can compare expected energies to minimal energies. 



1.1. Minimal Energies. 

First, we review results on CP"'^ = S^. We caution that different authors use 
slightly different normalizations of the radius of the sphere and the choice of Riesz 
kernels (l2| and the coefhcients of the asymptotic energies reflect these choices. 



In [KS jlBBPllRSZ] . the authors consider the energy ^ with 



Gsiz,w) 



z,w I 



{s y^O); lesp. Go{z,w) ^\og- (s = 0). 

z, w\ 



Here, [z, w] is the chordal distance between two points on S'^ C rather than the 
geodesic distance as in The are related by [z,w] = 2sin| on the sphere of 
radius | (i.e. area tt). In Si IWTl IW2| . it is proved that 



(5) 



s = 2 
s > 2, 



^V2is)N'^-C3{s)N^+i <£s{N)<^V2{s)N^^C4is)N^+i, s<2 



where all Cj{s) and ¥2(3) are positive constants independent of N. 

The asymptotic energy £o{N) on the standard was further studied in |BBP1 
IRSZj . Theorem 3.1 of [RSZ| asserts that if Cn is defined by 



(6) 



SoiN) = -] log(-)7V2 _ ^ logTV + CnN, 
4 e 4 



then for the standard S"^ one has 

TT 



1, nVs 

' - log = 

4 2 8^3 



> lim sup Cat > lim inf Cat > — log 

Af-i-oo Af-i-oo 4 



(1 



-a\h 



for certain explicit constants a, h. But the precise value of C^r is not known. 
In higher dimensional sphere S*™, results of [Bl IKS| show that 



(7) 



where Vm{s), jm and Cj{s) are positive constants independent of N. Many further 
references and results can be found in jBBGKS] . 

In the case of the Green's function G{z, w) of any compact Riemannian manifold, 
i.e. the kernel of on the orthogonal complement of the constant functions, N. 
Elkies proved the general lower bound 

^GreenW > -\^ogi^jN^ - ^NlogN - + oiN). 

Proofs can be found in [Hr) ILa[ lHiS| . In complex dimension one, the lower bound 
follows from Mahler's inequalities [M]. Although the right side is very precise, we 
are not aware of a comparable upper bound in the literature. 

1.2. Expected Energies of zeros of random polynomial systems. 

We now consider the energies of random configuration of points which are defined 
as the zeros of systems of random polynomials (or sections). Let {L,h) {M,Ld) 
be a positive holomorphic line bundle equipped with a Hermitian metric h — e~'^ 
{(j) is the Kahler potential. Let H'^{M, L^) be the space of the global holomorphic 
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sections of the line bundle L^; put (In '■= dim H'-\M , ) . For notation and 
background we refer to |GH1 ISZ2j . We normalize the volume form by 

(8, *:=f^v\ /n"=i. 



M 



We denote the zero set of a full system {s^ , . . . , s^} of m polynomials (or sec- 
tions) in dimension m by 

(9) = {(zi,...,z™)eM:sf(zi,...,z'") = 0, Vj' = 1, . . . , m}. 
The Riesz energy ([T]) of a configuration of zeros is thus given by 

(10) £:,(z,«,...,,„)= J2 Gsiz\zn 

Z^,Z^£Z iv N 

As in [SZT| [SZ2l IBSZlj and elsewhere, we endow the spaces H"{M,L'^) with 
Gaussian probability measures. They are induced by the inner product defined as 
follows: Let on denote the A^th power of the Hcrmitian metric h on L, and 
define 

(11) (si, S2) = / h^{suS2) , Si, S2 e H°{M, L^) . 

Jm ™' 

The Hermitian Gaussian measure on H'^{M,L^) is the complex Gaussian proba- 
bility measure /i^ induced on H'^{M, L^) by the inner product ( 11 1: 

(12) d/,A^(,)^_e-M^rfe, s^J^c.^f, 

where {Si ^ . . . , S^^} is an orthonormal basis for H^{M, L^). 

The Gaussian measure d^^ induces product measures on {M , L'^ ))'' = 
H^{M, L^) X • • • X H°{M, L^) {k times). For simplicity of notation we continue 
to denote them by d^f^ , leaving it to the context to make clear how may factors 
of H^{M,L^) are involved. Product measure corresponds to choosing < m in- 
dependent random sections, hence a Gaussian random system. We define E^nEs 
to be the expected (average) value of the energy of the zeros of Gaussian random 
sections chosen from the ensemble [H^ [M , ) , jif^) . In the discrete case, 
(13) 

Before stating our results in the discrete case, we recall the results of [Z) lABS] 
when m = 1. In the special case of the standard metric on 5^ and with the 
logarithm Riesz energy, it is proved in [Zj lABS] that 
(14) 



^ - i A^ log AT - f + oiN), radius = i 



~llogil)N^~lNlogN+l\og{^jN + oiN), radius =1. 

In our notation, 5*^ is the Kahler manifold (CP^, w^g), where ujps (or hps) is the 
Fubini-Study metric and the holomorphic line bundle L is hyperplane line bundle 
0(1); the area of {C¥^,ujps) equals 1, hence is the round S*^ of radius ^. The 
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expected energy of N randomly chosen points (i.e. from the bmomial point process) 
equals — | \og{^)N^ + ^ \og{^)N. Thus, the — iVlog term reflects the lower energy 
of well-separated points. 

To be precise, Q. Zhong obtained general asymptotic results for any Kahler 
metric on any Riemann surface, where the energy was defined in terms of the 
Green's function, 

(15) Go{z,w) = -^\ogrg{z,w) + Fg{z,w) 

where Fg is the Robin function such that J^^ Go(z, w)uj{z) = 0. In the case of the 
Green's energy, Zhong proved that 

(16) E^.£^ = _i7v log TV _ ^ _ iv Fg{z, z)io + o{N). 

where J^^ Fg{z, z)uj is the Robin constant. The N'^ term vanishes since the integra- 
tion of the Green's function is 0. In the special case of the standard metric on S"^, 
Armentano-Beltran-Shub |ABS| later showed that the expected energy is exactly 
given by the formulae in (14) without the o{N) term. 

The asymptotic results of [Zj lABSj show that the expected energy of zero point 
configurations is asymptotically the same as the minimum energy of any configu- 
ration up to order N (it appears to remain unknown what the coefficient of the N 
term is for the minimal energy). See also |ZZej for a correction and addendum to 
the calculations in 

1.3. Main Results in the discrete case. 

In this article, we consider random Riesz s-energies for zeros of polynomial sys- 
tems and their analogues on a general compact Kahler manifold of any dimension m 
and for a certain range of s. In the discrete case we obtain the following asymptotic 
expansion for the average Riesz energy. 

Theorem 1. Let (Af, w) be an m- dimensional compact Kahler manifold and let 
{L, h) — )■ {M, uj) he a positive Hermitian holomorphic line bundle. Then the expected 
Riesz energies defined in (13) admit the following asymptotics, 

• IfO<s< min{2m,4}, 

p-i 

E^^nEs =ai(/i, s)iV2" + ajih, s)N^"'-3 

+ c„,(s)7V™+4 +0(iV"+^(logAr)'"-i) 

where p = [m — |] + 1, where aj{h, s) are the coefficients in Lemma^ and 
where 

/•OD 

c„,(s) - 2m / [K„,^{r) - 1] r^^-^-^dr. 
Jo 

• In the case of logarithmic energy when s = 0, we have, 

E^nSo = ai(/i,0)A^2m ^ a2(/i,0)7V2m-2 ^ . . . ^ [aUh,0) - c^liV™ 

" h 

- iV™ log %/7V + 0(iV"-3 (logiV)™+i), 

where 

(17) c™ = 2m/ logrK„„,(r)-l]r2'"-idr. 
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Above, the coefficient c„i{s) involves the scahng hmit Kmm{T) of the pair corre- 
lation function of zeros in the discrete case in dimension m. As reviewed in |2.3[ it 
has the small distance asymptotics, 

(18) A^™™(r) = "^r4-2™ + 0„(r8-2-), as r -> . 

Hence the integral defining the coefficient Cm(s) only converges if s < 4. It is 
independent of the Hermitian metric h. 

The coefficients aj{h, s) are integrals of Gs{z,w) against curvature invariants of 
h and can be explicitly computed by the algorithm in Lemma [T] In particular, the 
first two are given by 



ai{h,s)^ Gsiz,w)dViz)dViw) 
(19) ^^'"*f 

a2(/i, s) = - 2 / Gs{z, w)A^S{z)dV{z)dV{w) 



J MxM 

where is the complex Laplace operator with respect to the Kahler form w, and 
S is the scalar curvature of w. 

We can give a more explicit evaluation in the special case of {C¥"^ ,Ldps) and 
{L,h) = {0{l),hps). As is well known, in a suitable frame, we can identify 
H°{CV"\0{N)) with the space of polynomials of m variables of degree N [GH] . 

Corollary 1. For the space H'^{C¥"^,0{N)) of polynomials of degree N equipped 
with the Gaussian Hermitian measure induced by hps, the Riesz energies of zeros 
of m independent random polynomials are given by 

• IfO<s< min{2m,4}, 

E^.£, = a^ihps, s)7V2m + c^{s)N"^+i + 0(iV"+"^(logiV)'"-f ) 

• For the case of logarithm random energies, we have, 

E^nSo = aUhFs,0)N^"' - iV" log VlV - c,„Ar™ + 0(Ar™-5(logiV)'"+i). 

^ h 

Thus, all of the terms aj{hps,s) vanish for j > 2. The proof follows directly 
from the proof of Lemma [l] and based on the fact that the Fubini-Study metric is 
balanced, i.e. the Szego kernel is constant on the diagonal. 

The asymptotics reflect the local behavior (repulsion/attraction) of the zero- 
point configurations and the long-range order (uniform distribution) as follows: 

• The leading coefficient ai{h,s) is the "Poisson term", i.e. the scaling limit 

of the expected energy of a binomial point process where one chooses N™ 

points at random from M, i.e. a configuration at random from •* ~ 

M X ■ ■ ■ X M /Spfm (the Af™-fold symmetric product) with respect to dV. 
^ 

See 92.41 for further discussion. 



The coefficients aj{h,s) are global lower order corrections due to the non- 
constant curvature of the Kahler metric. 

The terms Cm{s)N"^^^ and — iV™ log y/N ~ CmN'"^ reflect the local behavior 
(i.e. separation, repulsion or clumping) of the zeros. They are derived 
from the integration of Riesz functions against the current of two point 
correlations near the diagonal of M x M. The separation of the zeros 
depends on the dimension m. 
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• The overall asymptotics reflects a competition between the global and local 
behavior of zeros. The larger s is, the more significant is the separation of 
zeros. When s > 2, the local behavior of zero-point configurations domi- 
nates the global behavior; when s < 2, the global behavior dominates the 
local one. 

The asymptotic average energy of zero point configurations may be compared 
with the minimum energy ([7| , keeping in mind that the expected number of points 
in the configuration of zeros is N"^. Of course, S'" is not a Kahler manifold for n > 3, 
so no exact comparison can be made between the minimum and average energies. 
The closest Kahler analogue of a round sphere 5^™ is CP™ with its Fubini-Study 
metric of constant positive curvature. For Riesz energies with < s < min{2m, 4} 
the results of Corollary [l] may be compared with the result for the standard S*^™ 
that 

The exponent s/2 arises because we have A^™ points in real dimension 2m. The 
leading coefhcients of the minimum and average energies are both given by the 
Poissonian term and may be considered as the same. Hence the difference between 
minimum energy and mean energy of zero point configurations is first detected 
in the subleading terms of order N™^i, which as discussed above is due to the 
separation between points. The coefficient Cj{s) is not known explicitly, but it is 
significant to compare the signs of Cj{s) and c(s). It is known that Cj{s) > 0, 
so that the sub-leading term of the minimum energy is negative. Intuitively, we 
expect that Cm{s) > for s near 4 when m > 3 due to the clumping of the zeros 
as compared to the repulsion that must take place in an optimal configuration. 



The exact formula for Kmm(?') is given in [BSZ3 and recalled in (46 1. It is very 
complicated and represents the pair correlation function for zeros of random entire 
functions in the Bargmann-Fock model, i.e. entire analytic functions on C which 
are in i^(C, e~'^' ). Uramii") blows up to -l-oo for r — ^ if m > 3, supporting the 
intuition that Cm{s) should be positive. But it also has negative values in a certain 
interval, which represents a kind of repulsion relative to the Poisson spatial process 
for intermediate values of r, and that complicates the discussion. 
However we can justify the intuition for s near 4: 

Proposition 1. For any m > 3, there exists s,„ G [0, 4) so that for s > Sm, 

/>oo 

c,„(s) = 2m / [K„„,(r) - 1] r^^-^-^dr > 0. 



We do not try to find a good value of Sm, which probably requires a computer as- 
sisted calculation, but are content with a simple argument that nevertheless clarifies 
the structure of Kmmif)- 

1.4. The higher codimension case. The methods extend readily to systems 
of fc < m polynomials (or sections), whose zero set is a complex submanifold of 
codimension fc, so we also include this case. If A C M is a complex subvariety of 
codimension k < m ~ 1, the Riesz functions have no point mass on the diagonal, 
and we define the continuous Riesz energies as, 

(20) £s{A)= I I Gs{z,w)dVk{z)dVk{w) 

J J AxA 
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where dVi- = — rrr is the volume form on A. 

The zero locus Z^n ... of k holomorphic random sections is an analytic sub- 
variety of codimension k where k < m. We define, 



where 

(22) = f Gs{z,w)dVk{z)dVk{w) 

^1 >^k 1 ' k 

We then have, 

Theorem 2. i/ie same assumptions as in Theorem^ 

• If < s < 2(m — fc), w;e /laue, 



'^'^ \ (to — K)\ 



p-1 

ai{h, s) A^2fe ^ ^ ^f.^^. ^-^^ s)N^''~^ 

J=2 



rf,„(fc,s)7V2'=-'"+3 +0(7V2fc-m+V(logiV)™-«^ 



where p ~ [m — ^] + 1 and dmik, s) is the integral defined in ( 77). 
• If s — 0, we have, 

/ m-k \ 2 

E^nEo = [ai(/i, 0)A^2'- + W2ik)a2ih, 0)N^''-^ + ■■■ + w^{k)a„,{h, 0)iV2'=-™] 



where d{k) and e{k) are constants depending on k. 
where aj{h,s) are defined in Theorem^ and Wj{k) is the weight, in particular. 

Remark 1. Note that conditions of s in the two theorems are different. This is 
because the behavior of the pair correlation functions of zeros near the origin are 
different when k ^ m and k = ni respectively. 

1.5. Fekete points. Although we do not have any results on it, we should point 
out that the separation and equidistribution of points in dimensions to > 1 are 
often measured by the energies of Leja and Zaharjuta, 

(23) log AAr(zi, . . . ,zd„) = logdet(Sjv(2;i,2:j)) 

where _Bjv is the Bergman kernel for degree N polynomials (the kernel of the or- 
thogonal projection onto the degree N polynomials). Here, dN is the dimension of 
the space of polynomials of degree N . In dimension one, this is a Vandermonde 
type determinant whose logarithm is closely related to the Green's energy of the 
configuration. But in higher dimensions, the Green's energy and log det (i?Ar(zi, Zj)) 
are no longer known to be closely related and probabaly are not. Moreover, one is 



forced to choose different numbers of points for the two different energies: for ( 23 ) 
one chooses dN points, but in the Green's energy problem in the point case the 
number of zeros of the system is Ci(L^)'" > djsi- We refer to jBBWN) for recent 
result on Fekete points in higher dimensions. 
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2. Background 

In this section, we recall some definitions in Kahler geometry |GH| IZe| and some 
properties about zeros of random holomorphic sections of positive holomorphic line 
bundles pSZl. ■BSZ2. .BSZS^ . 

2.1. Kahler geometry. 

Let (M, w) be an m-dimensional polarized compact Kahler manifold. Then there 
exists a polarization of (M, w), i.e., there exists a positive holomorphic line bundle 
(L, h) with smooth Hermitian metric h such that the curvature form 

= -asiogiieLii^ , 

satisfies the following condition 



-- 2 -^'^ 

Here, cl is a local non- vanishing holomorphic section of L over an open set U C M, 
and lleill/i = /i(ei,eL)^/^ is the /i-norm of ,GHJ. 

We denote by H°{M,L^) the space of global holomorphic sections of = 
L(E)- ■ -^L. The metric h induces Hermitian metrics h'^ on given by ||s^^||/jiv = 



|s|j^. We give H'^{M,L^) the Hermitian inner product (11). 

We define the Bergman kernel as the orthogonal projection from the integral 
sections to the holomorphic sections; 

(24) Un{z, w) : L^{M, L^) ^ H°{M, L^) 

with respect to the inner product (•, 

Denote {5^ , • • • , 5'^} as the orthonormal basis of H^{M,L^) with respect to 
the inner product where — dim i?°(M, L^). Then along the diagonal, 

we have the following complete asymptotics jLl [Ze] , 



(25) HAr(z, z)=Y. Il^f = ^" + aMN"'-^ + a2(z)iV'"-i • • • 

where each coefficient aj{z) is a polynomial of the curvature and its covariant 
derivatives. In particular, ai{z) is the scalar curvature of w. 

Example: In the case of (CP™, wf^) with the line bundle (©(iV), hps) which is the 
A^th power of the hyperplane line bundle 0(1), the Bergman kernel is a constant 
along the diagonal [BSZ2j . 



(26) 
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2.2. Zeros of Random sections. 

For Gaussian random liolomorpiiic sections s^{z), ■ ■ ■ ,s^{z)e {{H'^{M, L^))'^, d/ij^) 
witli k < m, we denote by 

(27) Z,«^... {z e M : (z) = ■ ■ ■ = (z) = O} . 

It is a subvariety of codimcnsion k. The current of integration over the zero locus 
is defined by 

(28) - / 0(z), 0eI?"-'=-™-'=(M), 
which we also write as 

(29) =^^Z,«,...,,«A0(z), ^€V"^-'^^"^-\M) 
For the special point case of fc = m, 

(30) = 0(^), </'eC(Af). 

We now equip the vector space H^{M, L^) with the complex Gaussian probabil- 
ity measure ( 12 1 induced by h. This Gaussian is characterized by the property that 
the 2dN real variables 3?Cj, 3cj (j — 1, . . . , cZjv) are independent random variables 
with mean and variance |; i.e., 

(31) Ecj = 0, EcjCk = 0, EcjCk = Sjk ■ 

where E denotes the expectation 

We denote by Z^n the current of integration over the zeros locus of a single 
section. By the Poincare-Lelong formula. Then we have |SZlj . 

(32) E{Z,n) ^ -^ddlogllN + —uj 



It follows by (|25|) that 

(33) E{^Z,n)^- as N 



Define 



m 



(34) EZ^:^E{Z,.^...^,.), k< 

as the expectation of zeros locus of k random sections. Then we have |SZ1) . 

(35) = (EiZ,.))" = (^^dd\ogUN{z,z) + 
By ([25]) we get, 

(36) i5(z^„..._^„)^_^'= + 0(7v™-i) 

for N large enough. As a direct consequence, we know the expected number of 
discrete zeros is = (ci(L^))'" when k = m. 

Example: In the case 0{N) CP™ equipped with the Fubini-Study metric, 
H'^{CP"\0{N)) is identified with the space of polynomials of degree N. We equip 
H'^{C¥"\0{N)) with complex Gaussian ensemble induced by Fubini-Study met- 
ric. Since the Bergman kernel is constant (26), by Poincare-Lelong formula (32), 
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the distribution of zeros of Gaussian random polynomials satisfies the exact 
formula, 

(37) ^(i^V)-^ 

as N large enough. And the expected number of zeros is A^™. 

2.3. The pair correlation current. 

The pair correlation current is defined by, 

(38) K^rai^.w) := E ,,«(z) ® Z,«,... ,,«(«;)) e V^'''^\M x M) 

in sense that for any test form 0i(z) ® 02(w) G X>™-fc.™-fc(M) ® p™-fe,™-fe(M), 
we have the pairing, 

(39) {K^^{z,w),4>i{z)®^2{w))=E[(z^,,^...^,.,ct>i) , </>2) 
In the special case when k — m, the correlation measures take the form 

(40) K^.miz,w) = [A] A (ii;,^„(z)® 1) +^^(z,u;X®< 

where [A] denotes the current of the integration along the diagonal and E 
C°°[M X M) and E^^ is the expectation E{Z^n ... ^^n). 

The correlation currents of two points are smooth forms away from the 

diagonal in M x M and has no mass on the diagonal for A: < m; in the case of 
k = m, contains a delta-function along the diagonal |SZ1| . 

Iff < k < m — 1, the correlation current for is an current on AI x M 
given by 

(41) K^n = {-didid2d2QN + EZ,n EZ.Nf 

where Q^r is the pluri-bipotential defined in [SZlj . 
In the case of A; = m, 

(42) = {~didid2d2QN + EZ,« ® |mxm\a + diag, {EZ,«r 

where A is the diagonal in M x M and diag : M ^ M x M is the diagonal map 
diag{z) — [z, z). 

Furthermore, the pair correlation currents has the off-diagonal asymptotics [SZlj . 



(43) = {EZ,. ® EZ^^)" + 0{N-^), for rg{z,w) > 

where b > y/q + 2k + 3 > 3. 



2.3.1. Universal rescaling property. 

We now re-scale the correlation functions around a point zq G M in normal 
coordinates by a factor of In |BSZ2| , the authors prove that the pair correlation 
currents have the universal rescaling asymptotics. 



where the error term is times a smooth volume form. The universal rescaling 



property is independent of the manifold M and the line bundle L; depends 
only on the dimension m of the manifold and the codimension k of the zero set. 
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Furthermore, K'^{z,w) is a function depending on the distance rg{z,w) and thus 
we rewrite |BSZ3) . 

(45) K^{z,w) ^ Kkm{r), r=|z-w|, l<k<m 

The universal scahng hmit Kkm{T) is expHcitly computed in |BSZ21 IBSZ3) . In 
the case when k = m, one has the rather comphcated formula {v — ) 
(46) 

, , m(l-t>'"+^)(l-t>)+r^(2in+2)(i)'" + ^-i;)+r*[«'"+^+f™+({m+l}i)+l)()j"'-t))/(i)-l)] 
Kmm{r) = m(i_„)m+2 , 

For small values of r, we have 

(47) K^Ur) = I^r^-^™ + ©(r^-^™) , as r ^ . 

As mentioned above, the Kii(r) ~ as r — >■ 0, exhibiting repulsion, while K22(?') 
I when m — 2. For m > 3, K,mm{T) blows up as r — > 0. 
In the higher codimension case k ^ m, 

(48) Kfcm(r) - r-2'=, as r ^ 
In both cases, 

(49) Kfem(r) = 1 + 0(r*e"''^) , as r -J> oo . 

for any 1 < fc < m. Thus, the correlations are very short range on the length 
scale. The constant 1 would be the correlation function for the Poisson spatial 
process (i.e. independent points chosen at random from the normalized volume 
form). Below are computer graphics of Kmm{r) in the cases m — 2, 3, 4. 
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From the computer graphic, one sees that K33{r) — 1 and K44 — 1 each two zeros, 
both less than 4, and each is positive elsewhere. In Proposition [T] we show that 
kmm{f) — 1 > when m > 3 for > 2to + 3, but almost certainly it becomes 
positive at much smaller values of r. 

2.4. Configuration spaces, empirical measures and random points. Generic 
configurations of zeros C = {C^j • ■ • 1 C^"} i'^i the discrete (point) case are points in 
the punctured product 

(50) Mr, = {(C\ . ..,Cn^Mx---xM-.C^C' for p ^ qj/S,,, 



where = (ci(i^))'" is the number of zeros of the system and Sr„ is the sym- 
metric group on tn letters. As in the case m — 1, one has a map 

He : Mr, ^ Vm, fJ-C = :^Y1 
from configurations to probability measures on M. We also have the map ^ 



(5f , • • • , 4) e (H°(A./, L^))™ ^ e M, 



from polynomial systems to their zeros, and the composition gives the empirical 
measure of zeros map 

(sf ,sZ)^ f^c- 

The Gaussian product probability measure on (iJ'^(M, L^))™ pushes forward to 
a probability measure on Air, and the two-point correlation function (38) is the 



expected value of (g) with respect to this pushforward. Note that the image 
in Mr, of the zeros map ^ has positive codimension growing linearly in N, and 
therefore the induced measure on Mr, is very singular. 

By comparison, just choosing rjv random points on M from the normalized 
volume form Q corresponds to product measures on Mr, , known as the binomial 
point process. If we let — t- cx) and rescale as above, we get the Poisson point 
process on C" with intensity given by the Euclidean measure. Its correlation 
function is then the constant function 1. Hence Kmm ~ 1 repressents the deviation 
of the zero point correlations from the Poisson case. 

The Riesz energies of configurations ([T]) are functions on Mr, x Mr, which 
blow up on the 'large diagonals' where some pair — C'' (* 7^ j)- For polynomial 
systems, this large diagonal defines the discriminant variety 'Drn,N of full polynomial 
systems of degree N m m variables which have a double root. Thus £s defines a 
simple notion of 'inverse distance to 'Dm,N '• The large expected value for in 
dimensions m > 3 thus show that configurations of zeros are often close to Vm.N 
by comparison with random points. 



3. Proof of Theorem [T] 

Now we compute the Riesz energy of the discrete set Z^n ... which is the zero 
set of m random sections. In the proof, c is a constant, but may be different from 
line to line. First we need the following Lemma, 
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Lemma 1. We have complete asymptotics, 



(51) Jmxm 

= ai{h, s)7V2™ _|_ ^^(/^^ s)Ar2™-2 + ^^ih, s)N^"'-^ 

where each coefficient aj{h,s) is a constant depending on s and the curvature of h, 
in particular, the first two terms are given in uT 



Proof. From (35), E^^{z) which is the expectation of discrete zeros of m random 
sections is given by the exphcit formula, 

(52) e:^„,{z) = f-^aaiogn^(z,z) + -c.(z) 



First by (|25|), we have, 

1 



(53) lognAr(z,z) = mlogiV + ai(z)7V~^ + 



N- 



where ai = S{z) is the scalar curvature of uj and each aj{z) is a polynomial of cur- 
vature and its covariant derivative. Thus we get the following complete asymptotics 
oiEg^[z), 



(54) 



Y TT TT 



which implies that the integration (51) admits complete asymptotics. Each Ojih, s) 
is well defined since Gg G L^{M x M) for < s < 2m. 

We can in fact compute each coefficient from the asymptotic expansion of 
From the complete expansion of ( [54| , and the identity 

m{iddS) A uj""-^ = (A„S')a;™ 



we obtain ( 19 ) □ 
Now we are ready to prove Theorem [T] 

3.1. The case of < s < 2m. 

We first prove the case < s < 2m, the logarithm case when s — can be 
derived by modifications. By the definition of the discrete random Riesz energy, we 
have 

E^nSs E^N (^Gsiz,w), Z^N ... ,iv(z) (g) ... ^siv(w) - ^a(^)) 
= {Gs{z,w),KZm{z,w) - [A] A «„(z) ® 1)) 

G.(z, w) [kZ^{z, w) - [A] A ® 1)] 

MxM 

, Gs{z, w) [K^,^{z, w) - [A] A {E^,„,{z) ® l)] 

N 



, Gs{z,w)K^^{z,w) 

MJr,(z.w)>>^^^ 



By Lemma jlj we rewrite E^nEs as: 
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E^.Ss^ G.(z, w) [K^^„,iz, w) ~ [A] A {E^^„,iz) (g> l)] 

Jm ■/rg(z,tu)< ''^';^" 

, Gs{z,w)K^^iz,w)~ [ G,{z,w)EZ^{z)AEZM 

p-1 

lih, s)Ar2™ ^ J2 aj{h, s)N^"'-^ + 0(iV2™-P) 

Gs{z, w) [K^^iz, w) ~ [A] A (E^^^iz) ® l) - E^^{z) A 

p-1 

+ ai (/i, s) A^2m ^ s)iv2'"--'' + 0(iV2"-P) 

:=! + 11 + III 

where p is a positive integer which will be chosen later on. 
The following Lemma asserts that // could be ignored. 

Lemma 2. 

II (7V-«+*) 
where q is a positive integer which can be chosen large enough. 



N 



Proof. If rg{z,w) > , by the definition of Riesz functions, we have, 



\G,{z,w)\<c 



log TV 



as N large enough. If we choose fc = m in the estimate (43), 

K^^{z,w)-E^^{z)AE^^^iw) = 0{N-^') on r,iz,w)>^'^^ 



where b > ^/q + 2k + 3. If we combine these two estimate, the term // reads, 
G,iz,w) [K^U^,w) - eZ^{z) a <„(ii;)] = 0{N-'i+i) 



which completes the proof. □ 

Now we turn to estimate / to finish the proof of Theorem [l] 
Lemma 3. 

/ = c„(s)A^™+5 + O (^iV'"+ ^ (log 0<s<2m,s<4, 
where c(s) is a constant given explicitly in the proof. 
Proof. We change variable 

u 

w = z ^ = 
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Then by the rescaUng property of (44), we have, 

- [A] A «,„(z) ® 1) =KZ^{z,z+^) - [A] A {EZ^(z) ® 1) 

(55) 



where Kmm(|u|) is defined in (|45| with k = m. 

'uj{zy 



By the estimate of ( 54 1 , 
(56) E^,,^{z)AE^Jw)^N'"^ 



By choosing the normal coordinate, if Lp is the Kahler potential, then we have 



27V Vi 



Since u — v']V(w — z), we have 



Therefore, 
(57) 



dy,dyjip{w) = dudu^iz + ^) = dudj-^ + 0{N ^ ) 

V N V 



Thus we rewrite (56 1 near the diagonal when \u\g < fo-^log N as, 



(58) E^,^{z) A £;^„(z + ^) - (^) 



27r 



Recall the definition of /, 



M Jrg{z,w)< ' 



Gs{z,w) [kZ^{z,w) - [A] A {E^^iz) l) - <„(z) A <„(ii;)] 



Now we apply identities (|55|), (|58| to rewrite 

I = N 

(59) 



G,(z,2 + -=) 

M Ja<\u\g<b^\ogN VJV 



|u|)-l + 0(iV-5) 



- A 



2tt 



dd\u\ 



The error term 0{N 2) comes from estimate (551 and (58). 
Thus we further rewrite ( 59 ) as, 

/ 



(60) 



Gs{z, Z + -—) [K„„(|li|) - 1] 



27r 



dd\u\ 



Denote I = I + Rg- The estimate of the error term Rg is as follows. 



Rs = 0(iV™-3) 



Gs{z,z+ —= 

M J\u\g<b^\ogN vN 



27r 



(99|u| 
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If we choose the polar coordinate, we have 



M Jo<\u\g<b^/^^^ 
< c [ r^H^"'-'d\u\ < ciVi(logiV)"-i 

In the last step, we apply the estimate of the geodesic distance function rg{z,w) 
near the origin, 

(61) rg{z,z+^) = ^+0{N-i) 



in the normal coordinate. Thus, 

(62) i?,, = O (^iV"+^(logiV)"-4) if0<s<2m, 

The estimate of the integration / is as follows, first we rewrite, 

JmJg<\u\<ooH ^tt/ \2t: J 

+ lower order terms 



by applying the expansion of the geodesic distance function (61). We use polar 
coordinates to compute the leading term. 



^K™(H)-i]0"A(^aa| 



J M J0<|«|<oo F, 
'»+*^2m-l [ 1 r _ 11 l„|2m-l 



=A^"+^ / ^ K™(kl) - 1] \ur''d\u\ 

(53) ^ Jo<|«|<oo Fl 

= ( 2m / [Kram{r) " 1] r^^-^-^dr ) iV™+i 

\ "'o<|m|<oo y 

:=c„,(s)7V™+f 

where V2m-i is the volume of (2m — l)-dimensional unit sphere. 

The integrand in c(s) is always integrable on the tail a < |u| < cx) as a large 
enough, this is easy to see because of asymptotic property Km„i(r) = 1 + 0(r'*e~'' ) 
as r — > oo. But we need some conditions to ensure c(s) to be well defined near the 
origin. Since we know Kmrnif) n^iir^"^™ as r — )■ 0, it's easy to check that s < 4 
is a necessary condition. 

Moreover, we have the easy estimate, 

Bs = ©(TV^+i-i) 

Thus the growth of the error term Rg will dominate the growth of Bg, and 

I = As + Rs^ Cm{s)N"'+i + O (7V™+'^(log7V)'""5) 
which completes the proof. □ 
We now finish the proof of the main result Theorem [T] for < s < 2m. 
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Proof. For the case < s < 2m and s < 4, recall the decomposition of E^nEs, 

p-i 

E^nEs = ai(/i, s)iV2™ + aj{h, s)N^"'~^ + 0{N'^"'^p) + I + II 

3=2 

Now we apply the estimates in Lemma [2] and Lemma |3j we have, 

p-i 

E^nEs ^ai{h, s)7V2" + J2 aj{h, s)7V2m-j _^ QiN^^'-P) 

+ c™(s)A^™+3 +o(iV™+^(logiV)™-3) +0(iV-«+3) 

The smallest integer p such that 2m — (p— 1) > to+| will ensure the above 
expansion is well defined, i.e, p = [m — |] + 1. The last term is negligible for q large 
enough. Thus we rewrite, 

p-i 

E^nEs =ai(/i, s)iV2™ + J2 aj{h, s)N^"'^^ 

+ c„^{s)N"'+i +0 (Ar™+'^(logAr)'"-f^ 
which completes the proof. □ 

3.2. The case of s = 0. 

Now we sketch the proof for the case of s = in Theorem [T] 
First, as before, we rewrite, 

E^.Eo = - [ [ \ogr,{zM {Kl^{zM - [A] A ®\)~ E^^^z) A EI^{w)\ 

/ , logrg(z, w) [K^m{z, w) - EZra{z) A 

p-1 

+ ai {h, 0)iV2™ + J2 aj [h, 0)N^"'-^ + C»(7V2™-p) 

:= I + 11 + III 
Follow the steps in Lemma [2j 

// = 0{N~'^ log A^) for q large enough 



since IGqI is bounded by logA^ if rg{z,w) > ''^^'^ . Thus the term II can be 
ignored in the final expansion. 

The statement in Theorem [l] about logarithm Riesz energies follows directly by 
the following lemma. 

Lemma 4. 

(64) I = -N"" log Vn - c^N"" + ©(iV^-s (log iV)™+i) 

where c is a constant given by [P, 
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Proof. Follow the proof of (60) in Lemma [s] we rewrite / = 7 + i?o as follows, 



(65) 



u 

Go{z, Z + ——) [Kmr„(|u|) - 1] 



(- 



a(— +i?o 



where 



Using the geodesic distance estimate (61|, wc have, 



/o<|u|g<6VlogJV 



log ( % ) \u\^^-^d\u\ 



< c(log A^) 



m+1 



Thus, 



i?o - O (7V™-5(log7V)™+i) if s = 0. 
We further rewrite I = Aq + Bq by applying the expansion of the geodesic distance 



function (61) again, 
Ao=-N' 

J M ^0<|m|<oo 

(66) =iV"logV]V 
- I 2to / 



log 



u 
N 



[K.nun{\u\) - 1] 



|Kmrn(|u|) - 1] ( —dd\u\ 



/0<|m|<oo 



0<|«|<oo 

logr[«;„,„(r)-l]r2'"-idr ) TV' 



Now claim: 
(67) 



/ 

Jo< 



[Kmm(|w|) - 1] ( 7;-dd\u\ 
\ Zn 



-1 



J0<\u\<oo 

We calculate the integration directly by the rescaling property: 



M ^""^ JO<|u|<c 

1 c r 



[Kmm{\u\) - 1] [^9d\u\' 



— lim , , 

N^oo TV™ Jm Jr,(z,w)<'- 



{kZ^{z,w) - [A] A £,^^(1.) - eZ^{z) a i?^„(z«)} 



: lim 



JM Jrg{z,w)>- 

The claim follows by the following two calculations: By the definition of correlation 
currents, the first term reads, 

N^co iV™ y^/iAf 

= lim —{N^"" - N"" - N^"") = -1 

N^oc, iV™ ^ ^ 
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The second term reads, 



lim — — 

lim 0(iV-«-") ==0 



where we apphed estimate (43) for q large enough. Thus 



(69) 



Ao = -iV'" log - CnN 



where €,„ is given by (17). The integral above is also well defined by the asymptotic 



properties of Hmmi^)- And it's also easy to check, 



Bq = 0{N"'-') 

The growth of Bq is less than the growth of Rq, thus the error terms Bq can be 
absorbed in the error term Rq. Thus in the case of s = 0, 



1 = Ao + Ro = -A^™log c^TV" + 0(iV'"-5(log7V)™+^) 



which completes the proof. 



□ 



4. Proof of Proposition [T] 

The proof is almost trivial but we include it to clarify the structure of Kmm(r) — 1 
(46). The main point is that Cm(s) is a meromorphic function with poles at s = 
4, 2m. When m > 2, the only pole for s G [0, 4] is at s = 4. When m > 3 the 
residue is positive, while for m = 2 it is negative. When m = 1 it also has a pole 
at s = 2 S [0,4] with a negative residue. The Proposition just uses these facts to 
draw conclusions on positivity near enough to the poles. 

The results are illustrated in the following computer graphics in the cases m = 
2, 3, 4. The graphs for to > 4 are rather similar to that for m = 3, 4. The integral 
in (71 ) or more precisely {Kmm ~ l)r^"^~^~^dr is evaluated numerically over the 
interval [0, 10000] so that positivity of the numerical integration implies positivity 
of the true integral. As the figures show, Cm(s) < for small s, in which case the 
random zeros have smaller Riesz s-energy on average than general random point. 
When 771 = 2 it appears always to be negative, so that random zeros have less 
energy than random points on average. When to > 3, Cm(s) > for s somewhat 
larger than s — 1 and the rise is very steep as s — >■ 4. As to increases, the clumping 
of zeros increases and hence the Riesz energies increase. 
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12 3 




m=3,4 

Here are the details: We observe that the numerator of Kmm consists of three 
terms 









= m(l - 




(70) 


< 




= r2(2m + 2)(^;™^ 












f ({m + l}w + - v)/{v - 1 



Since v < 1 only the second is negative. Furthermore the second and third terms 

2 

are divisible hy v = and therefore the first term dominates for large r. We 
claim that there exists ro(m) so that Kmm(r) — 1 is positive for r > ro(m). We do 
not attempt to find the best ro(m) but rather look for a simple one. 

We first drop the terms u™ or higher powers from the numerator , simplifying 
the numerator to 
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with X = r^. We would like to find x so that No{m,x) > m{l — e i.e. so 

that 

o.(m+l) _^ 1, e^^ (2m + 2) _„ _„ ,^ -TMn-i-? 

x^[- -e ^ -xe ^-e ^>(l-e ^)'"+-^ _ i. 

m m 1 — e ^ m 

In fact, the left side is easily seen to be positive if x > (2m + 3). Returning to our 
original numerator, we want to show that NQ{m,x) > m ii x > (2m + 3). Given 
the above observations, it suffices to show that the missing terms — mw'"+^ (1 — v) + 
r^(2m + 2)w™+^ > and this also holds when a; > 3. Consequently, Kmm{r) — 1 > 
if a; > (2m + 3). Thus, we can take ro{m) = yj (2m + 3). Hence, to prove that 
Cm(s) > if suffices to show that 



/.V2m+3 

(71) / [K„„,(r)~l]r2™-i-^dr >0 
Jo 

for s > s,„. Indeed, the contribution for small r has the form, 

i-e 777 + 1 f^^^ ,2m-s 

(72) / [^™™(r)-l]r2"-i-dr=^- +0,„(e«-^), 

Jq 4 4 — s Im — s 

and only the first term has a pole as s — > 4^ when m > 3. When m = 2 the integral 
of —1 also blows up at s = 4 and the coefficient of is —1/4, so the integral is 
negative for s near 4; when m = 1 the integral of —1 term tends to — oo at s = 2. 
Hence in proving positivity we assume m > 3. On the other hand, for m > 3 



the integral of the left side of ( 71 ) over the interval [e, ^/2m + 3] is continuous in 
s € [0, 4] and of course has a uniform upper bound. Hence ( |7l| ) holds for s < 4 
sufficiently close to 4. 

5. Proof of Theorem [2] 

The proof of Theorem [2] is similar to Theorem [T] For simplicity, we only prove 
the case when s ^ 0, the result for logarithm energies can be derived by few 
modifications. We now sketch the proof as follows. 

Proof. Suppose Z^n ... is the zero locus oi k < m holomorphic random sections, 
i.e., analytic subvariety of positive dimension m — k, hy the definition of the con- 



tinuous random Riesz energy \21\ and ( 22 1 , we have, 
Gs{z, w)dVk{z)dVk{w) A K^^(z, w) 

MxM 

Gs{z,w)dVk{z)dVk{w) A i?L(z) A E^M 

MxM 

+ f Gs{z, w)dVk(z)dVk{w) A [K^^iz, w) - E^^{z) A E^M] 

J MxM 

=:I + II 

m-k 

where dVh — — rxx- 

{m—k)\ 
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Following the step in Lemma [ij E^^{z) which is the expectation of zero locus 
ZgN ... ^jv admits the following expansions, 



(73) 



--N 



-TT/ TT 



fc-1 



iddS{z) 



Thus, it's easy to get, 



. , Mh,s)N'^ +W2{k)a2{h,s)N"''' +w:i{k)a:i{h,s)N'''-' + ■■■] 

[m — k)\ J 

where all aj are the same as in Lemma [T] and Wj is the weight of a^, in particular, 

h 

^ m 

To the second term //, we need to decompose it to be. 



M Jrg(z.w)<- 

III + IV 



M J Tg {Z.'W)'> - 



Recall the estimate (43), 



Xf„(z, w) - E'^^(z) A Eg^{w) = OiN-") on r,(z, w) > 



this is true for any k. Thus following the argument in Lemma [2] 

IV = 0(iV-«+3) 

for q large enough. Thus IV can be ignored in the final expansion of E nEg. Thus 
we will finish the proof if we can get the estimate of ///. 
We change variables, 

u 

w = z A = 

Then by the rescaling property of (44), we have, 

(74) i^L(z, ^ + = N'^^uM) (^)' A {^dd\u\''f + 0{N^-^ 



Following the step in proving the estimate ( 58 ) , we have 

k 



(75) Et,{z)^Eg,X^+—)=N^' 



A I —dd\u\^ 
Zn 



Thus we further rewrite /// as. 



/// = N' 



2fe-r 



i — k 



J- — Ml) / / Gs(^,2+ ^) [Kfe„(|w|) - 1] 



A I ^dd\u\^ 



where 

Rg = 0(Ar2fe-m-i) 



g,(z,z + ^)(^)"a f aah 
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Thus by the same argument in Lemma [3] we have, 

(76) i?3 O (7V2fe-™+^(iog7V)™-5) if0<s<2m, 



Applying the asymptotics of geodesic distance function (611 near the origin, we 

rewrite, 

(77) 

/// ^N'^-rn+i f f ^ [.U\U\) - 1] (-)" A f f 

\{m-ky.J 7m Jo<|u|<oo |u| ^tt/ \2tt 

+ Rs + lower order terms 



\{m-k)\) vr™ Jo<|u|<oo 
+ -Rs + lower order terms 
= : dm(s, fc)iV^'''""+5 + i?g + lower order terms 

by polar coordinate. 

The integrand in (i(s, fc) is always integrable on the tail a < |u| < cx) as a large 
enough because of asymptotic property Kkm{i") = 1 + 0{r^e^^ ) as r — > oo. To 
ensure dm{s, k) to be well defined near the origin, recall that when k ^ m, 

Kkm{r) ^ r^^'' as r 0, 
hence, s < 2(m — fc) is a necessary condition. 

□ 

Now we can turn to prove Theorem [2) 
Proof. For the case s 7^ 0, s < 2(m — fc), recall the decomposition of E n£s, 



2 



{m-k)\, 

+ d„,{s,k)N'^^-"'+i +0 (iV2''^-™+'^(log7V)™-5) +0(7V-«+4) 

The smallest integer p such that 2fc — (p — 1) > 2fc — m + | will ensure the above 
expansion is well defined, i.e, p — [to — |] + 1 which completes the proof. If we 
further choose q large enough, we obtain 

Ea^£s = 7 TTT [ai{h, s)N^'' + W2{k)a2{h, s)N^''-^ + ■■■] 



□ 



(to — fc)! ^ 

+ rf„,(s,fc)A^2fc-m+f +0 (^7v2fc-™+^(iogAr)™-f j 



The same method is also applied to logarithm energies when s = 0, but we omit 
the proof here. 
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